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Abstract 

We derive the main physical galaxy properties: mass, halo radius, phase space density and veloc- 
ity dispersion from a semiclassical gravitational approach in which fermionic WDM is treated 
quantum mechanically. They turn out to be fully compatible with observations. The Pauli Prin- 
ciple implies for the fermionic DM phase-space density Q{r) = p(r)/<x 3 (r) the quantum bound 
Q(r) < K m 4 /h 3 , where m is the DM particle mass, cr(r) is the DM velocity dispersion and K 
is a pure number of order one which we estimate. Cusped profiles from A^-body galaxy simu- 
lations produce a divergent Q(r) at r = violating this quantum bound. The combination of 
this quantum bound with the behaviour of Q(r) from simulations, the virial theorem and galaxy 
observational data on Q implies lower bounds on the halo radius and a minimal distance r m! „ 
from the centre at which classical galaxy dynamics for DM fermions breaks down. For WDM, 
r m i„ turns to be in the parsec scale. For cold dark matter (CDM), r m ,„ is between dozens of kilo- 
meters and a few meters, astronomically compatible with zero. For hot dark matter (HDM), r„„„ 
is from the kpc to the Mpc. In summary, this quantum bound rules out the presence of galaxy 
cusps for fermionic WDM, in agreement with astronomical observations, which show that the 
DM halos are cored. We show that compact dwarf galaxies are natural quantum macroscopic ob- 
jects supported against gravity by the fermionic WDM quantum pressure (quantum degenerate 
fermions) with a minimal galaxy mass and minimal velocity dispersion. Quantum mechanical 
calculations which fulfil the Pauli principle become necessary to compute galaxy structures at 
kpc scales and below. Classical A^-body simulations are not valid at scales below r mm . We apply 
the Thomas-Fermi semiclassical approach to fermionic WDM galaxies, we resolve it numeri- 
cally and find the physical galaxy magnitudes: mass, halo radius, phase-space density, velocity 
dispersion, fully consistent with observations especially for compact dwarf galaxies. Namely, 
fermionic WDM treated quantum mechanically, as it must be, reproduces the observed galaxy 
DM cores and their sizes. The lightest known dwarf galaxy (Willman I) implies a lower bound 
for the WDM particle mass m > 0.96 keV. These results and the observed galaxies with halo 
radius > 30 pc and halo mass > 4 x 10 5 M Q provide further indication that the WDM particle 
mass m is approximately in the range 1-2 keV. 
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1. Introduction and summary of results 



Dark matter (DM) is the main component of galaxies, especially of dwarf galaxies which are 
almost exclusively formed by DM. It thus appears that studying galaxy properties is an excellent 
way to disentangle the nature of DM. 

Cold DM (CDM) produces an overabundance of substructures below the ~ 50 kpc till very 
small scales ~ 0.005 pc which constitutes, as is well known, one of the most serious drawbacks 
for CDM. On the contrary, warm DM (WDM), that is, DM particles with mass in the keV scale, 
produces DM structures in the observed range of scales > 50 kpc. In WDM structure formation, 
substructures below the free-streaming scale ~ 50 kpc are not formed contrary to the case of 
CDM. This conclusion for WDM based on the linear theory is confirmed by A^-body simulations 
1 11,0 SHE 53, Sill]- For scales larger than 50 kpc, WDM yields the same results than 
CDM and agrees with all the observations. 

Astronomical observations show that the DM galaxy density profiles are cored till scales be- 
low the kpc 1 71 , 61 , 5^, 15 , 24 , 65 ] . On the other hand, A^-body CDM simulations exhibit cusped 
density profiles with a typical l/r behaviour near the galaxycenter r — 0. WDM simulations 
exhibit cusps or small cores smaller than the observed cores 

B0I 



6J,|64J] 



Numerical calculations on the spherically symmetric Vlasov-Poisson equation based on the 
Larson moment expansion 134 13311 . as well as on the exact dynamics of the associated N-body 
system, have confirmed these findings 11711 . 



For fermionic DM the Pauli principle states that the phase-space distribution function for 
spin-i particles f(r, p) must be smaller than two 



f(r,p)<2 



(1) 



Since the matter density p{r) is obtained from the phase-space distribution through 



(2 n hy 



m = DM particle mass 



(2) 



this implies a bound on the phase-space density Q(r) = p(r)/cr 3 (r) where cr{r) is the DM velocity 
dispersion, which takes the form (see sec. |2]i 



— Off) < K . 
nr 



(3) 



Here A" is a pure number of order one that we estimate in appendix Appendix A and display in 
Table Q] In the classical physics limit H — > 0, the right hand side goes to infinity and the bound 
on Q(r) disappears. 

The quantum bound eq.(Q]i in the cosmoiogicai context has been previously considered in 
refs. lf59l EM 0] to derive lower bounds on neutral lepton masses. In the present paper, we 
use this bound in a new way incorporating the behaviour of the phase-space density Q(r) |Q1 
I27L I37I l62L I45II and observational data to derive lower bounds on the halo radius. These bounds 
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are of semiclassical gravitational nature and motivate us to use the Thomas-Fermi approach for 
galaxies, namely, treating gravitation classically and fermionic WDM quantum mechanically. 

We combine the quantum bound eq.([3]l with the results from classical A^-body simulations 
which indicate a simple power law behaviour for Q(r) ll, 27, 37, 62, 451 



Q(r) = Qh [y\ (4) 

where j3 1.9-2, is the halo radius and Qi, stands for the mean phase space density in the 
halo. Classical WDM A^-body simulations yield results similar to the behaviour eq.(H]leven for 
r < r%, when properly using the velocity dispersion of the keV-scale WDM particles 11381.15 ill . 

Classically, Q(r) grows unbounded for r — > violating the quantum bound eq.©. Therefore, 
for small enough scales classical physics in galaxies breaks down. In other words, classical N- 
body simulations and hydrodynamical simulations are not applicable near the galaxy centers. 

We show that the classical result eq.© only holds for distances far enough from the galaxy 
center, namely for 

r > r min = — \ — r h ■ (5) 
m- V K 

For a DM particle mass in the keV scale and the observed values of r/, and Qh, it turns out that 
r m t„ can be between 0.1 and 1 pc. 

Eq.© provides the minimal distance from the center where quantum mechanical effects are 
essential and rules out the presence of galaxy cusps for fermionic WDM. Moreover, quantum 
mechanical effects extend well beyond r = r m ,„ because smoothing up the classical cusp which 
sucks matter towards the origin has an effect on the whole galaxy halo. 

Taking into account the quantum nature of fermionic DM introduces a physical length scale 
of quantum origin. As a consequence, density profiles become regular (cored) at the origin. For 
bosonic DM, the bound eq.dOb does not apply and the formation of cusps is allowed. However, 
the observed galaxy density profiles are cored 10 EDS [II Hill. 



The quantum bound eq.(TT3l applies for any kind of fermionic dark matter with the value 
of r,„i„ determined by the DM particle mass. For CDM where 1 GeV < m < 100 GeV, r m ,-„ is 
between dozens of kilometers and a few meters, astronomically compatible with zero. Therefore, 
classical approaches as A^-body simulations and the classical Boltzmann-Vlasov equation fully 
apply for CDM and unavoidably produce cusps. 

HDM where 1 eV < m < 10 eV, suppresses all structures scales below the Mpc scale due to 
its long free-streaming length and has been ruled out years ago. Anyway, we find that classical 
aproaches to fermionic HDM as A^-body simulations are not valid for scales below an r,„,„ which 
turns to be between the kpc and the Mpc, depending on the galaxy type. 

A direct way to see whether a system of particles has a classical or quantum nature is to 
compare the particle de Broglie wavelength A^b with the interparticle distance d. We do that in 
sec. 13. H and express their ratio as 
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The observed values of Qh from Table [2]yields 



2xl(T 3 <<R i-^r-X < 1.4. (6) 
\kev/ 

The larger value of % is for ultracompact dwarfs and the smaller value of % is for big spirals. 
The values of % around unity clearly imply (and solely from observations) that compact dwarf 
galaxies are natural macroscopic quantum objects for WDM. 

WDM fermions always provide at least a pressure of quantum nature. When this quantum 
pressure is balanced with the gravitation pressure, we find values for the total mass M ~ 10 6 M Q , 
the radius R ~ 30 pc and the velocity dispersion <x ~ 2 km/s consistent with compact dwarf 
galaxies (see Table [2). These results back the idea that dwarf spheroidal galaxies are supported 
by the fermionic WDM quantum pressure. 

We then treat the self-gravitating fermionic DM in the Thomas-Fermi approximation. In this 
approach, the DM chemical potential p(r) — po - m <p(r), where po is a constant and <p{r) the 
gravitational potential, obeys the Poisson equation 

d 2 p 2 dp 4nGm 2 f°° 2 

+ --r = TTT~ P d P f ( e (p) ~ P( r )> (') 



dr 2 r dr n 2 h? 



where G is Newton's gravitational constant, p is the DM particle momentum, e(p) = p 2 /(2m) is 
the DM particle kinetic energy and f(E) is the energy distribution function, 

This is a semiclassical gravitational approach to determine selfconsistently the gravitational 
potential of the fermionic WDM given its distribution function /. 

In order to have bounded DM mass densities according to eq.© we impose the boundary 
condition at the origin: p'(0) = 0. 

The distribution function f(E) is determined by the DM evolution since decoupling. Such 
evolution must take into account the quantum character that emerges when the distribution func- 
tion approaches the quantum upper bound eq.Q. Such quantum dynamical calculation is beyond 
the scope of the present paper. We modelize here the distribution function f(E) by the equilib- 
rium Fermi-Dirac distributions and by out of equilibrium distributions for sterile neutrinos (see 



the appendix Appendix A i 



We get a one parameter family of solutions of eqs.((7} parametrized by the value of the chemi- 
cal potential at the origin /i(0). We then express the chemical potential at the origin in terms of the 
phase-space density at the origin Q(0). Large positive values of //(0) correspond to the quantum 
degenerate fermions limit, while large negative values of p(0) yield the dilute (classical) limit. 
We show that the Thomas-Fermi equation implies the local equation of state 

P(r)= i v 2 (r)p(r) , 

and the hydrostatic equilibrium equation 

dP dd> 
-r;+p(r)-f=0. 
dr dr 

This local equation of state generalizes the local perfect fluid equation for r-dependent velocity 
v(r). 
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In the quantum degenerate fermions limit, the halo radius, the velocity dispersion and the 
galaxy mass take their minimum values. These minimum values are similar to the estimates for 
degenerate fermions given in sec. 13.21 

The masses of compact dwarf spheroidal galaxies dominated by DM must be larger than 
this minimum mass M„„„. The lightest known galaxy of this kind is Willman I (see Table [2]). 
Imposing M„„„ < Mwuiman i — 3.9 10 5 M G provides a lower bound for the WDM particle mass: 

m > 0.96 keV . (8) 

The numerical resolution of eqs.(|7]i for the whole range of the chemical potential at the origin 
fi(0) yields the physical galaxy magnitudes: mass, halo radius, phase-space density and velocity 
dispersion fully compatible with observations especially for compact dwarf galaxies as can be 
seen from figs. [T]and|2]and Table|2] 

Approaching the classical diluted limit yields larger and larger halo radii, galaxy masses and 
velocity dispersions. Their maximum values are limited by the initial conditions provided by the 
primordial power spectrum which determines the sizes and masses of the galaxies formed. 

The phase space density decreases from its maximum value for the compact dwarf galaxies 
corresponding to the degenerate fermions limit till its smallest value for large galaxies (spirals 
and ellipticals) corresponding to the classical dilute regime. 

The theoretical values for r/,, M and v(0) from the Thomas-Fermi approach vary very little 
with the specific form of the phase-space distribution function f(E) as function of the energy 
(Fermi-Dirac distribution or out of equilibrium sterile neutrino distribution). 

Comparison of the theoretically derived galaxy masses with the galaxy data in fig. [^indicates 
a WDM particle mass m approximately in the range 1-2 keV. For larger masses, m » 1 keV an 
overabundance of small galaxies (small scale structures) without observable counterpart appears. 

In summary, the theoretical Thomas-Fermi results are fully consistent with all the observa- 
tions especially for dwarf compact galaxies as can be seen from figs. Q~|and[2] It is highly re- 
markably that in the context of fermionic WDM the simple static quantum description provided 
by Thomas-Fermi is able to reproduce such broad variety of galaxies. 

These results indicate that fermionic WDM treated quantum mechanically (even approxi- 
mately) is fully consistent with the observed galaxy properties including the DM core sizes. 

Therefore, the effect of including baryons is expected to be a correction to the pure WDM 
results presented in this paper, consistent with the fact that dark matter is in average six times 
more abundant than baryons. 

In section|2]we describe the quantum bound for the phase-space density of fermionic DM and 
we derive the relevant galaxy scales where DM classical physics in galaxies breaks down. In sec- 
tion III we show how the quantum fermionic WDM pressure balances the gravitational pressure 
in galaxies. In section IV the Thomas-Fermi approach for self-gravitating fermionic WDM is 
presented and applied to galaxies, showing that fermionic WDM treated quantum mechanically 
is able to reproduce the observed DM properties of galaxies including the DM cores and their 
sizes. We derive in Appendix A the numerical values of the constant K in the quantum bound for 
several momentum distributions. 

We use units where the speed of light is taken c = 1 along this paper. 
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2. Quantum bounds for fermionic DM in galaxies from the Pauli principle 

For fermionic DM the Pauli principle tells us that the number of spin states in a phase- 
space cell volume (2n ft) 3 cannot be larger than 2 for spin-| particles. Namely, the phase space 
distribution function f(r, p) must satisfy 

f(r,p)<2. (9) 
The DM number density can be expressed as 



_ f(r,p) = nf_ _ 
(2 Tift) 3 2 ft 3 



= d P T^rir = TTT <r(f) f(fi K , (10) 



where f(r) is the maximal p-average of the phase space distribution over a volume m 3 cr 3 (r), 
m is the mass of the DM particle, cr(r) is the DM velocity dispersion, <r 2 (r) =< v 2 (r) > /3 and 
K is a pure number of order one. In appendix Appendix A we evaluate K for distributions of 
cosmological relevance and display it in Table [J 
We find from eq.(TT0t 

2 ft 3 n(f) 

/W=t^-^t (11) 
Inserting eq. dTTb into eq.© which also applies to f(f) yields, 

n{f) ^ Kn? (12) 



(T 3 {?) ~ ft 3 

As a consequence, we find the following bound for the phase-space density: 

where p(r) = m n{f) is the matter density. Therefore, the phase space density Q{r) can never 
take values larger than the right hand side of eq.dTjt. This is an absolute quantum upper bound 
which is due to quantum physics, namely the Pauli principle. 

In the classical physics limit ft — > the right hand side of eq.([T3ll tends to infinity and the 
bound disappears. 

Since the squared velocity dispersion cannot be larger than the speed of light we have a lower 
bound for the phase-space density 

Q> Qmin = 3 V3p . 

The upper bound eq.([T3li on Q implies, at given density p, a lower bound v„„„ on the velocity 

Qmax = 3 V3 = K ^ . 

mm 

and therefore 

i 

v > v min = ft V3 (jr-A 3 (14) 
6 



Distribution Function 


1 


K 


0.60364 

rrz 

VK 


Fermi-Dirac and DW 


3 £(3) = 3.60617 


1.89858 


0.43809 


v-MSM 


0.045 n f(5) = 0.14659 


0.077178 


2.17285 


Maxwell-Boltzmann 


y|= 1-2533 


0.65985 


0.74311 



Table 1: The constant K denned by eq.JB} and 0.60364/ We see that the characteristic size is larger for sterile 
neutrinos decoupling out of equilibrium. 



N-body simulations yj, |27J, [371 |62j, |45|] as well as the resolution of Larson's equations IU7I1 



point to a cuspy phase-space density behaviour 



Q(r) = Q h y-\ (15) 

where f3 1.9-2, ry, is the halo radius and Qi, stands for the mean or characteristic phase space 
density in the halo. Classical WDM /V-body simulations yield results similar to the behaviour 
eq.© even for r < r/, when properly using the velocity dispersion of the keV-scale WDM 



particles 0381,1511 



We see that this Q(r) derived within classical physics tends to infinity for r — > violating 
the Pauli principle bound eq.([T3l. Therefore, classical physics in galaxies breaks down near the 
galaxy center. 

Setting /? = 2 for simplicity, we find by combining eqs.(fT3l and (fT3T > 

r > r min = — \\ — r h . (16) 
m- \ K 

That is, r > r„,j„ given by eq.(fT6b sets the domain of validity of classical physics in DM dominated 
galaxies. 

The quantum bound eq.([T3l) rules out the presence of galaxy cusps for r < r,„,„ and implies 
the existence of galaxy cores with size larger than r„„„. 

The values of r/, and Qi, vary several orders of magnitude according to the type of the galaxy 
(see Tabled. Qh is larger for dwarf spheroids than for spiral galaxies. From Table|2] 

10-4 <^LVSa< 1.3 (17) 
(keV) 2 
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The larger value corresponds to ultra compact dwarf spheroidals and the smaller one to spiral 
galaxies. We finally have from eqs. ( TToT l and (fTTI i, 



l(r 



keV\ 2 13_ /keV\ 



(18) 



r/, goes from ~ 10 pc for ultra compact dwarf spheroidals till the 10 kpc range for spiral galaxies. 

The right hand larger value in eq.([T8~t corresponds to the smaller galaxies (ultra compact 
dwarfs) with r/, in the 10 pc range. The left hand smaller value in eq.(fT8l corresponds to the 
larger galaxies (spirals) with r;, in the 10 kpc range. Multiplying both sides of eq.(fT8l by r/, we 
obtain for a DM particle mass in the keV scale: 

0.1 pc < r min < 1 pc , (19) 

because r/, is larger for the more dilute galaxies where Qi, is smaller and vice versa. 

Notice that eq.dToli provides the minimal size of the region where quantum mechanical ef- 
fects are strong. Quantum mechanical effects will influence a region well beyond the lower 
bound r,„i„. Namely, removing the classical cusp which sucks matter towards the origin has an 
effect on the whole galaxy halo. 

The quantum bound eq.© in the cosmological context has been previously considered in 



ie q 

refs. |59[ |25|, |26|] to derive lower bounds on neutral lepton masses. In the present pap er, we 
incorporate to this bound the behaviour of the phase-space density Q(r) (Hill, 37, 62[ 45] to 
derive lower bounds to the halo radius eq.([T6V These bounds are of semiclassical gravitational 
nature. 



Assuming, as usual, that virialization holds within the radius r/,, namely 

2 _ GMj, 

where Mi, is the mass within the radius r/,, 



(20) 



M h = \n r \p h , Ga=3, (21) 
3 h 

and G is Newton's gravitational constant, we obtain from eqs. (fToT ), d20l and (fSTJ, 

- — I 3 

Qh ~ 4n \M h {Gr h f ' 



(22) 



and 



h y 3 r h V 3 h ImpTf I 3 r h V 



2 m 2 \ VG/ \ M h I 2 yf^K \ m ) \G M h 

where mpi = y/h/G is Planck's mass. More explicitly 

r min 0.58793 I r h 1O 6 M \^ /keV " 



P c \P C M h 



(23) 



(24) 



V_ J ill LlA V 


n 


<T 




o(0)/ M ° 


M h 


pc 


km 


(keV) 2 


(pc) 3 


10 6 M Q 


Willman 1 


33 


4 


0.85 


6.3 


0.39 


Segue 1 


38 


4 


1.3 


10 


0.59 


Leo IV 


151 


3.3 


0.2 


.19 


1.14 


Canis Venatici II 


97 


4.6 


0.2 


0.49 


1.43 


Coma-Berenices 


100 


4.6 


0.42 


2.09 


1.97 


Leo II 


233 


6.6 


0.093 


0.34 


7.25 


Leo T 


152 


7.8 


0.12 


0.79 


7.4 


Hercules 


305 


5.1 


0.078 


0.1 


7.5 


Carina 


334 


6.4 


0.075 


0.15 


9.6 


Ursa Major I 


416 


7.6 


0.066 


0.25 


12.6 


Draco 


291 


10.1 


0.06 


0.5 


21 


Leo I 


388 


9 


0.048 


0.22 


22 


Sculptor 


375 


9 


0.05 


0.25 


22.5 


Bootes I 


322 


9 


0.058 


0.38 


24 


Canis Venatici I 


750 


7.6 


0.037 


0.08 


27.7 


Sextans 


1019 


7.1 


0.021 


0.02 


35 


Ursa Minor 


588 


11.5 


0.028 


0.16 


56 


Fornax 


944 


10.7 


0.016 


0.053 


74 


NGC 185 


355 


31 


0.033 


4.09 


293 


NGC 855 


837 


58 


0.01 


2.64 


2480 


Small Spiral 


4800 


40.7 


0.0018 


0.029 


5100 


NGC 4478 


1490 


147 


0.003 


3.7 


1.96 x 10 4 


Medium Spiral 


1.73 x 10 4 


76.2 


3.7 x 10~ 4 


0.0076 


6.4 x 10 4 


NGC 731 


4850 


163 


9.27 x 10~ 4 


0.47 


8.52 x 10 4 


NGC 3853 


4110 


198 


8.8 x 10~ 4 


0.77 


8.54 x 10 4 


NGC 499 


6070 


274 


5.9 x 10~ 4 


0.91 


3.27 x 10 5 


Large Spiral 


5.18 x 10 4 


125 


0.96 x 10~ 4 


2.3 x 10~ J 


5.2 x 10 5 



Table 2: Observed values for r h , a, V3fc, p(P) and M h for galaxies from refs. (65|, H^, [T^, H3 153|, HE S H3. - 
The phase space density is lar ger for smaller galaxies both in mass and size while the surface density /.( = p(0) r;, is 
approximately constant [22, 30, 57, 19]. Notice that the phase space density is obtained from the stars velocity dispersion 
which is expected to be smaller than the DM velocity dispersion. Therefore, the reported Qi, are in fact upper bounds to 
the true values {53j. 
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where we used that 



2/ V^5? (k pc )| Ml Gi (keV) 

The values of the constant 0.60364/ are displayed in table Q] while those of r^/pc and 
10~ 6 M h /M B are displayed in table |2] 

It must be stressed that r,„,„ is the minimal value for the core radius. The core radius can 
be well above this lower bound which corresponds to maximally packed fermions around the 
center of the galaxy. Hence, one expects for diluted objects as galaxies core radii much larger 
than the lower bound r„„„. Even in atoms the phase-space density turns to be significantly below 
the quantum bound eq.([T3l for thermal fermions rf32ll . Moreover, our derivation of 77, within the 
Thomas-Fermi approach presented in sec. |4]shows that r/, » r,„,„ for dwarf galaxies. 

The quantum bound eq.dT3t applies for any kind of fermionic dark matter. However, the 
value of r m i„ strongly depends on the DM particle mass. 

For CDM where 1 GeV < m < 100 GeV, we see that r„„„ can be from eq.(l24li between 
dozens of kilometers and a few meters. That is, r,„,„ for CDM is astronomically compatible with 
zero. Therefore, classical aproaches as A^-body simulations and the classical Boltzmann-Vlasov 
equation fully apply for CDM and unavoidably produce cusps. 

For HDM where 1 eV < m < 10 eV, we find r m ,-„ from eq.d24Ti between the kpc and the Mpc, 
depending on the galaxy type. HDM suppresses all structures scales below the Mpc scale due to 
its long free-streaming length and has been ruled out years ago. Anyway, it is important to learn 
that classical aproaches to HDM as A^-body simulations are not valid for scales below r mm which 
can be from the kpc to the Mpc scale according to the type of galaxy and the value of m in the 
eV scale. 

For bosonic DM, the bound eq.(fT3l does not apply and the formation of cusps is allowed. 
However, the observed galaxy density profiles are cored Jill III EI Hill (see also Q). 

In all cases cusps of fermionic DM in the galaxy density profile are artifacts produced by 
classical physics computations irrespective of the nature of dark matter (HDM, WDM, CDM) 
and of whether the computations are numerical or analytical. 

That is, quantum physics, namely the Pauli principle, rule out galaxy cusps for fermionic 
dark matter. 

We have so far ignored baryons in our analysis of the galaxies . This is f ully justified for dwarf 
spheroidal galaxies which are composed today of 99.99% of dark matter 



large galaxies the baryon fraction can reach values up to 1 - 3 % J48 , 4a " 43ll . We have also 
ignored supermassive central black holes which appear in large spiral galaxies but have not been 
observed in dwarf galaxies. As it is known, the formation of supermassive central black holes 
is correlated to the formation of the galaxy itself but this whole issue is beyond the scope of the 
present paper. Anyway, it must be noticed that the central black hole mass is at most ~ 10" 3 of 
the mass of the bulge. Baryons can be important in the galaxy formation and evolution, a subject 
that is outside the scope of our paper. Fermionic WDM by itself produces galaxies and structures 
in agreement with observations. Therefore, the effect of including baryons is expected to be a 
correction to the pure WDM results presented in this paper, consistent with the fact that dark 
matter is in average six times more abundant than baryons. 
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2.1. Classical upper bounds from the DM phase-space density evolution 

The quantum upper bound on the phase-space density eq.(IT3T> is unrelated to the classical 
upper bound derived from the self-gravity classical dynamics as already noticed in ref. ll59ll . 
Namely, the fact that the coarse-grained phase-space density when averaged also in r-space over 
a galaxy halo, can only decrease by classical collisionless phase mixing or classical self-gravity 
dynamics 1281114 , 

SSSEMl- This averaged phase-space density Q is thus bounded by its 



primordial value Q pr j m . 

This upper bound follows from the classical self-gravity dynamics evolution contrary to 
eq.(fT3b which applies to fermionic DM and has a quantum origin on the Pauli principle. In 
addition, this classical upper bound only applies to the phase-space density averaged over a large 
volume and on time intervals: Q < Q pr im- 

At specific points like the galaxy center the phase-space density is classically unbounded due 
to the appearence of cusps and certainly much larger than the primordial phase-space density. 



The primordial DM phase-space density which is space-independent can be written as II18I1 



V27 gj 4 / 2 2 

Qpr "" = 2n^ m TT (25) 
'4 

where g, is the number of internal degrees of freedom of the DM particle (g, = 4 for Dirac 
fermions), I2 and I4 are the dimensionless momenta of the primordial DM distribution function. 
Using eq.(f25ll for Q pr im, its self-gravity dynamical evolution and observational data for the phase- 
space density today (as in Tabled indicated that the DM particle mass is in the keV scale fljll . 

Inserting eq.(l25b into the quantum bound eq.(IT3b yields 

V27 l\ < K n 1 1\ . 

It is easy to check that this inequality is fulfilled for the distribution functions considered in the 



appendix Appendix A as it must be. 

That is, the primordial value Q prml of the phase-space density safely fulfills the quantum 
bound eq.(fT3b. 



3. Gravitation and the WDM pressure: the main drivers of galaxy formation 

As is well known, galaxy formation as all structure formation in the Universe is driven by 
gravitational physics. 

The main notion regarding structure formation at a given scale r is the Jeans' length Aj(f) at 
this scale. The Jeans' length is the product of the speed of sound or the velocity of the particles 
v times the free fall time 

Ur) = v(r) J-pArx (26) 

where G is Newton's gravitational constant. 

Length scales larger than Aj(r) are gravitationally unstable. Therefore, Aj(r) indicates the 
scale of the largest objects that can be formed at the distance r. The mass density p(r) decreases 
from its maximum value at the center of the galaxy r = till asymptotic typical values 200 times 
the average density in the Universe at ^200 ~ Rviriai- 
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At a scale r structures larger than Aj(r) can not form because they are unstable while struc- 
tures smaller than A j(r) can grow reaching a size ~ r. Therefore, 



r ~ Aj{r) . 



During the linear regime of the cosmological evolution the phase space density is small 
enough and the Pauli principle is always fulfilled. Notice that the phase space measure d 3 p d?r 
is invariant under the universe expansion and therefore the universe expansion by itself does not 
affect the fullfiling of the Pauli principle. 

Only when structure formation increases the particle density and particle velocity, the phase 
space distribution function f(r, p) may increase to approach its quantum upper bound 2 for spin 
one-half fermions. 

In summary, starting from the linear cosmological fluctuations at WDM decoupling, we can 
classically evolve the phase space distribution function f(r, p) as long as / remains well below 
its quantum upper bound. As soon as classical physics breaks in some region, even small, we 
have to use a quantum mechanical treatment not only in the region where a classical / would 
violate the quantum upper bound but in the whole region around. This is not an easy problem 
and may be solved using quantum Monte Carlo methods, time dependent Hartree-Fock methods 
and quantum Boltzmann equations evolution 



Clearly, one can use the classical evolution in most of the space where / stays well below the 
quantum bound and match it with the quantum evolution around the galaxy centers. 



The Jeans' mass is given by 
M 

In terms of the phase-space density, A j and Mj read 



<27) 



^y = 3i J-=—= , Mj = 



4 n~- 



GvQ ' 3 J \G 3 e 

Then, from eq.([3]l and eq.© with ft — 2 we have the bound 



Aj(r) - \nrh Vf (28) 

y K cr(r) V G mr 

More explicitly, we get a lower bound of the product of the Jeans' length and the square root of 
the velocity dispersion 



r— 0.339 /keVr 



showing that neither the Jeans' length nor the velocity dispersion can vanish for r — > 0. 
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3.1. Dwarf galaxies as WDM quantum macroscopic objects 

To determine whether a system of particles has a classical or quantum nature we should 
compare the particle de Broglie wavelength with the interparticle distance. 

The de Broglie wavelength of DM particles in a galaxy can be estimated as 

Ada = — , (29) 
m cr 

while the average interparticle distance d can be estimated as 

'-£)' ' 

where pi, is the average density in the galaxy core. By using pi, = o 3 Q/, and eqs.d29l-(l30ll. we 
can express the ratio 



as. 



* 5 d 



K = H(%) 3 . (31) 
Using now the observed values of Q/, from Table [2]yields 

4 

2x 1(T 3 <<R — < 1.4 (32) 
VkeV/ 

The larger value of % is for ultracompact dwarfs and the smaller value of % is for big spirals. 

The ratio % around unity clearly implies a macroscopic quantum object. 
Notice that here as well as in the bound eq.([T3l fi 3 Qjm A measures how quantum or classical 
is the system (the galaxy). 

We conclude solely from observations that compact dwarf galaxies are natural macroscopic 
quantum objects for WDM. 

3.2. Dwarf Galaxies supported by WDM fe anionic quantum pressure 

For an order-of-magnitude estimate, let us consider a halo of mass M and radius R of 
fermionic matter. It can be fermionic DM or baryons. Each fermion can be considered inside a 
cell of size Ax ~ 1/m and therefore has a momentum 

h * i 
p fin 3 . 

Ax 

The associated quantum pressure P q (flux of the momentum) has the value 

Pq = no~p~ho~n'i = — n s . (33) 
m 



where cr is the mean velocity given by 



p h i 
cr — — = — n 3 . 
m m 
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The system will be in dynamical equilibrium if this quantum pressure is balanced by the gravita- 
tional pressure 

G M 2 1 

Pa — gravitational force/area = — — x (34) 

5 R 2 4nR 2 

We estimate the number density as 

M 



3 

and we use that p — m cr to obtain from eq.(l33l 



4 ttR 3 m 



3 M 
R 5 \4n m 



P q = —. It— I ■ (35) 



Equating P q = Pq from eqs.(t34b-(l35l> yields the following relations between the size R and the 
velocity cr with the mass M of the system: 

3i h 2 /1O 6 M \^ /keV\^ 

R= 5 r = 10.6 ...pc — - (36) 

(4 np G m~ M' \ M ) \ m ) 

(AnV G 4 : km / m \l / M V 

cr = \ — — m~ M~ = 22.9 . . . — — . (37) 

h s IkeV/ \1O 5 M / ' 

Notice that the values of M, R and cr are consistent with dwarf galaxies. Namely, for M of the 
order 10 6 M (typical mass value for dwarf spheroidal galaxies), R and cr have the correct order 
of magnitude for dwarf spheroidal galaxies for a WDM particles mass in the ke V scale (see Table 
0. 

These results back the idea that dwarf spheroidal galaxies are supported by the fermionic 
WDM quantum pressure eq.( l35b . 

It is useful to express the above quantities in terms of the density p, as follows 



M=^ ^ = 0.7073... 10* M Jp^(^) 



«- ^=^=31.05... pc (-*)' (SX 



^2 P^ onn m-ll / pC 3 \ 3 /keV 



CT = 



f> = i— = 4.399 10" 1 ' ^ — , (40) 

q m l p C 3 \ p M J \ m j 

R and M are typical semiclassical gravitational quantities involving both G and h. The particle 
velocity cr and the pressure P q are of purely quantum mechanical origin. 

The radius R and mass M are the semiclassical Jeans' length and Jeans' mass as can be seen 
by inserting the velocity from eq.(|39l in eqs.(l26li and (l27T i with the result 

2n /2tt\ 3 
h= — R , Mj= — I 
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The associated phase-space density results 

4 

m 

Q = Jf 

In this case where v has purely quantum origin [eq.((39}], the de Broglie wavelength is equal 
to the interparticle distance Ajb = d. Thus the ratio % is unity which is the extreme (fermion 
degenerate) quantum case. 

In the presence of squared angular momentum L 2 we have to add the centrifugal pressure 

1 L 2 



4ttR 2 M R 3 

in the equilibration equation P q + Pi — Pq- The relation between the radius, mass and velocity 
eq.d36b takes now the form 

G M 3 (4 g m- M~ 
A simple estimation of the total angular momentum of the whole halo goes as follows 

L 2 ~ - M 2 R 2 3 a 2 
2 

where the \ factor comes from averaging the sin 2 of the angle between the momentum p and the 
particle position f, and the factor 3 comes from the relation v 2 = 3 a 2 . We thus obtain 



R= 10.6... pel " I I— V +3.48...pc 



M \ ml M 



R 



10 10 pc 



2 



(42) 



We see that the angular momentum contribution increases the size R. However, for dwarf galax- 
ies, we see that R and <x have the same order of magnitude for L > and for L = 0. 

4. Galaxy properties from quantum fermionic WDM in the Thomas-Fermi approach 

DM particles are nonrelativistic during structure formation and their chemical potential is 
given by 

/"to = Mo ~ m <Kr) (43) 
where m is the mass of the DM particle, jjQ is a constant and <p{r) is the gravitational potential. 

We consider for simplicity the spherical symmetric case where the Poisson equation for the 
gravitational potential takes the form 

^ + 1 = -4nGmp{r) (44) 
dr z r ar 

where G is Newton's gravitational constant and p(r) is the DM mass density. 
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Since the DM mass density is bounded by the Pauli principle as analyzed in sec. [2] p(r) must 
be bounded at the origin and therefore we must impose as boundary condition at the origin: 

^(0) = . (45) 
ar 

We can write the DM mass density as the integral over the momentum of the DM distribution 
function 

PW = Ala P 2 dpf[e(p)-^r)] (46) 
2 n z nr Jo 

where e(p) - p 2 /(2 m) is the particle kinetic energy, f(E) is the energy distribution function and 
g is the number of internal degrees of freedom of the DM particle, g — 1 for Majorana fermions 
and g — 2 for Dirac fermions. 

Eqs.(l44b and ( l46l i provide a system of ordinary nonlinear differential equations that deter- 
mine the chemical potential jj(r) and constitutes the Thomas-Fermi semi-classical approxima- 
tion. Fermionic DM in the Thomas-Fermi approximation has been previously considered in ref. 

guanas. 

The Thomas-Fermi equations provide a semi-classical approximation in which the energy 
distribution function f(E) is given. 

The boundary condition for the chemical potential is usually set to zero at a given radius 



R which is the appropriate condition for the Coulomb interaction in atoms IB2I1 and for self- 
gravitating degenerate fermions 113 ill . In the selfgravitating degenerate case the mass density is 
proportional to the power 3/2 of the chemical potential and therefore both quantities vanish at 
the same point. In the general selfgravitating case for non-degenerate fermions where eq.(l46t 
applies, the mass density stays nonzero beyond the point where the chemical potential changes 
sign and becomes negative. 

We integrate the Thomas-Fermi nonlinear differential equations (t44b -(f46b from r — till the 
boundary r — R — R200 ~ R V ir defined as the radius where the mass density equals 200 times the 
mean DM density. 

In order to determine the distribution function f(E) one should follow the precise DM evolu- 
tion since decoupling. Such evolution must take into account the quantum character that emerges 
when the distribution function approaches the quantum upper bound eq.©. This quantum dy- 
namical calculation is beyond the scope of the present paper. 

We modelize here the distribution function f(E) by equilibrium Fermi-Dirac functions and 
by out of equilibrium Dodelson-Widrow and v-MSM distribution functions (see the appendix 



Appendix Ai 



Let us define dimensionless variables v(£), £ , £0 and the dimensionless distribution function 
^as 



E 



(47) 



where Eq is the characteristic comoving energy of the DM particles at decoupling. The charac- 
teristic length Lq emerges from the dynamical equations (I44l -(|46]> and is given by 



VG(2m) 2 \E 
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For g — 2 the mass density takes the form 



= ^ hf) , I3(v) = / dy *¥(y 2 - v) 



(49) 



where we use the integration variable y = p/ y/2 m Eo and the Poisson equation (f44b-d46T> become 

d 2 v 2 dv 
+ ? dg 



— + ~ — = -3 flvtf)) , p(fo) = 200 Pdm , AO) = , (50) 



The distribution functions *F are given in the Appendix Appendix A for the thermal and out of 
equilibrium cases. 



In order to integrate eq.(l50l we have to specify v(0). v(0) is determined by the value of the 
phase space density at the origin Q(r = 0). That is, to find a given galaxy we must give Q(r = 0) 
and the boundary radius R in dimensionless variables £o. In eq.d66b below we find the relation 
between v(0) and Q(r = 0) which permits to compute v(0) from a given value of Q(r = 0). 

4.1. Main physical galaxy properties: mass, velocity dispersion, density and pressure 

We derive in this subsection the expressions for the main physical properties of the galaxies 
(mass, velocity dispersion, density and pressure) in the Thomas-Fermi semi-classical approxi- 
mation. 

The average velocity of the particles is space-dependent and follows from the average mo- 
mentum given by 

2 . . 1 f~p 4 dpf[e(p)-n(r)] 2E 2 

v (r) = — — — = a (v(£)) , (51) 

m j P 2 dp f[e(p) - n(r)] « 

where 



a(v) = 



J Q yUy^-v) 

A 75-1 : • (52) 

y 2 dy V(y 2 - v) 



The mass enclosed in the sphere of radius R follows by integrating the mass density given by 
eq.® 

M = 4* drp{r) = {^j £f d{ /3(v(0) , (53) 

where we used eqs.(l47b and ( l48l l. The integral over £ can be performed with the help of eq.dSOt. 

f # P(v(t)) = ~\& v '(&) ■ (54) 

Notice that v'(^o) is always negative since the function /3(v) is positive definite. 
We can write from eqs.(f53T>-(l54l> the mass M as 
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The pressure at the point r, P(r) is given by an analogous integral [Q] 

in 1 n? \ m } Jo 3 7T- \ m / 

From eqs.d49l and (l5TT i we can write the density and the average square velocity as 

p(r) = p(0) — — , v(r) = v(0) , 

po 

Here a = «(v(0)) and = /?(v(0)). 

Moreover, from eqs.(l49"V (IBTt and (l56l we derive as local equation of state: 

p (r) = l v 2 (r)p(r) . (57) 

This local equation of state generalizes the local perfect fluid equation of state for r-dependent 
velocity v(r). As we see below, the perfect fluid equation of state is recovered both in the classical 
dilute limit and in the quantum degenerate limit. 

It is very instructive to compute the derivative of the pressure eq.(l56l> with respect to r = Lq 
Upon integrating by parts and using eqs.(l43ll. (l47l i and (l49l we find 

dP dd> 

_+p(r)-f=0. (58) 
dr dr 

This shows that the quantum Thomas-Fermi equation implies the hydrostatic equilibrium equa- 
tion. 

By analogy with the Burkert density profile, the halo radius ry, = Lo & in this theoretical 
calculation can be defined in dimensionless variables as, 

Pfa) = K&) = I 
p(0) A) 4 ■ 

Furthermore, from eqs. (l49l and (l5"Tl the phase space density Q(r) is given by 

V27 4 #0 • ~ m V27 4 A) /r 
= TTT m ^77T and 2(0) = "~ 2~tt '« -r , (59) 

which turns to be independent of Eq. 

From eqs.d47li. d48l and (IB3T l. R 3 M turns to be independent of Eq too: 

, (3 nf h 6 s | , | 

We expressed above the main physical galaxy magnitudes Lo, Af, v(r), P(r) and p(r) in 
terms of the DM characteristic energy £o and the potential v(£), solution of eq.d50l). Because Eq 
is not directly observed we will eliminate Eq in terms of p(r = 0) from eq.j49l. [We may also 
choose other point r + 0]. We obtain 

(60) 



m J rrr j8o 
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By eliminating Lo and using eqs.d49li-(f52t. we have for the pressure eq 



P(r) = — h 2 m 4 
3 



P(0) 



a 2 (v) /3{v) and P(0) = y h 2 m 4 



\p(oy 


s 

3 


f \ 


m 4 







(61) 



Inserting eq.® in eqs.(|47]l, (US), (ED and {55]l yields 



3 n 

8G m \ 



P(0) 



v(r) = ft Q'(v) 



^ 2 P(0) 



More explicitly, we can express Lo, M and v(r) as 





PqM q ' 




p(0) pc 3 



3 ft 
8G keV - 



M 



V^ft 3 2 , , p(0) 
(62) 



7r pc 



= 22.47 pc , 



(63) 



Mo s 



Vgjri ft 3 
2§ keV 4 




1.425 10 3 M R 



v(r) = v 



keV 



«(v) 



vo s ^ 3 



^f/keV« 
pc 3 



= 4.268 ^ 

s 



The semiclassical galaxy magnitudes Lo, M, v(0) and P(0) emerging from the Thomas-Fermi 
equations generalize the corresponding expressions eq.(l38l-(l39]l derived just equating the grav- 
itational and WDM quantum pressures and describing the quantum degenerate fermions limit. 
Eqs.(l63l-(l64li cover the full range of physical situations from the degenerate fermions till the 
dilute classical limit as we discuss in the next subsection. 

The local Jeans' length and Jeans' mass follow by inserting eqs.d49l and (IBTT l for p{r) and 
v(r), respectively, into eqs.(l26ii and d27l i 



MO 



Lo 



a{0 



8tt 3 
V27 



M 



(65) 



We see that Aj and Mj differ from Lo = R/€o and M by ^-dependent factors of order one. 
Therefore, the galaxy length scale Lo and galaxy mass M emerging from the Thomas-Fermi 
equation < T50b are a measure of the Jeans' length Aj and the Jeans' mass Mj, respectively. 

From eqs.d49ll. ( 1521 and using Q(0) = 3 V3 p(0)/v 3 (0), eq.d64l) can be rewritten as 



a _ {ry 4 dyny 2 -y(0)]Y _ 2.145. 
4 {f~y 2 dyV[y 2 -v(0)]f h 



keV 4 



G(0) 



I m \i 

= 1.0772 

\keV/ 



M 



2(0) P c 3 (if) 



(64) 



/ m \~3 
\keV7 



(66) 

Eq.(l66b shows that v(0) is determined by the phase density at the origin Q(0). In the next sub- 
section we solve the Thomas-Fermi eqs.(l50ll in the whole range of the chemical potential at the 
origin v(0). 
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4.2. Physical galaxy properties from the resolution of the Thomas-Fermi equation 

Large positive values of the chemical potential at the origin v(0) » 1 correspond to the 
degenerate fermions limit which is the extreme quantum case and oppositely, v(0) <*: - 1 gives 
the diluted limit. The diluted limit is the classical limit and in this case the Thomas-Fermi 
equations (144-b - (f46b become the equations for a selfgravitating Boltzmann gas. 

The bounds on the phase space density follows from eqs.d49ll, (l52l and d59l ), 

, 2(0) 5 V5 > VO(0) 

< ft 3 ^ T-V = 0.37760 ... and < m , < 0.61449 ... for -oo < v(0) < oo 

(67) 

The largest value for the phase space density corresponds to the degenerate fermions limit while 
the smallest values appear in the classical dilute limit. 

In the quantum degenerate fermions limit, the halo radius, the velocity dispersion and the 
galaxy mass take their minimum values: 



24.516 pc | 



m 



P(0) ^ 



,(0) = 4.768 



(¥)* 


p(0) pc 3 







(68) 



These minimum values are similar to the estimates for degenerate fermions eqs.(l38)-(139K as it 
must be. 

The masses of compact dwarf spheroidal galaxies dominated by DM must be larger than the 
minimum mass M m - m eq.d68ll. The lightest known galaxy of this kind is Willman I (see Table|2|. 
Imposing M„„„ < Mwuiman i = 3.9 10 5 M gives a lower bound for the WDM particle mass: 

m > 0.96 keV . (69) 

Approaching the classical diluted limit yields larger and larger halo radii, galaxy masses and 
velocity dispersions. Their maximum values are limited by the initial conditions provided by the 
primordial power spectrum which determines the sizes and masses of the galaxies formed. 

The phase space density decreases from its maximum value for the compact dwarf galaxies 
corresponding to the degenerate fermions limit till its smallest value for large galaxies (spirals 
and ellipticals) corresponding to the classical dilute regime. 

Thus, the whole range of values of the chemical potential at the origin v(0) , given by the 
boundary condition, from the extreme quantum (degenerate) limit v(0) » 1 to the classical 
(Boltzmann) dilute regime v(0) <K -1 yield all masses, sizes, phase space densities and velocities 
of galaxies from the ultra compact dwarfs till the larger spirals and elliptical in agreement with 
the observations (see Table|2]i. 

In the degenerate limit the equilibrium FD thermal case and the out of equilibrium case give 
identical results, as expected. 
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In fig. Q]we plot from eq.d59ll the dimensionless quantity 




(70) 



In fig. |2j we plot the dimensionless product 



M Mr 




1.425 10 : 




■)| 



(71) 



M yjp(0) pc 



where M is the galaxy mass and we used eqs.(l63ll-(l64l. In both figures we plot in the abscissa 
the product 



where r/, is the halo radius. The phase-space density Q(0) and the galaxy mass M are obtained by 
solving the Thomas-Fermi eqs.dSOt for thermal (FD) fermions and for out of equilibrium sterile 
neutrinos with the distribution eq. (IA.6b . 

Notice that the theoretical curves in the right-hand side of eqs.(l70l-(l72l are independent of 
the value m of the DM particle and do not change with m. 

In fig. Q]we have superimposed the observed values of y/Qii/m 2 for m = 1 keV, m — 2 
keV and m - 5 keV (see Table Notice that the observed values of Qi, from the stars' ve- 
locity dispersion are in fact upper bounds for the DM Q/,. This may explain why the theoretical 
Thomas-Fermi curves in fig. [T]appears below the observational data. Notice that the error bars 
of the observational data are not reported here but they are at least about 10 - 20%. 

In fig. |2]we have superimposed the observed values of (M/M Q ) y/M Q /[p(0) pc 3 ] (m/keV) 4 
for m = 1 keV, m = 2 keV and m = 5 keV (see Tabled. 

We see from fig. [T]that increasing the DM particle value just pushes down and to the right the 
observed values ffl y/Qk/m 2 . In fig. |2]we see that increasing the DM particle value just pushes 
up and to the right the observed values (M/M Q ) y/M Q /[p(0) pc 3 ] (m/keV) 4 . 

As noticed above from eqs.(l70Ti-(l72t. the theoretical curves are independent of the value m 
of the DM particle. Hence, for growing m > keV the left part of the theoretical curves will have 
no observed galaxy counterpart. Namely, increasing m » keV shows an overabundance of small 
galaxies (small scale structures) without observable counterpart. This is a further indication that 
the WDM particle mass is approximately in the range 1-2 keV. 

The galaxy velocity dispersions from eq.(l64l turn to be fully consistent with the galaxy ob- 
servations in Table |2] 

We see in figs. [T]and|2]that fermions at equilibrium and out of equilibrium give very similar 
values for Q{0) and M. The theoretical values for r/,, M and v(0) vary very little with the 
distribution function *P. This is similar to the WDM linear power spectrum ||20ll , where changing 
the distribution function can be balanced by changing the mass m within the keV scale. 

Notice that the values obtained for the halo radius r/, in fig. Q] are much larger than the 
quantum lower bound r„„„ eq.dT9b, as expected. 

In summary, the theoretical Thomas-Fermi results are fully consistent with all the observa- 
tions especially for dwarf compact galaxies as can be seen from Table [2] This result gives an 




,3 



(72) 
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additional support to the idea put forward in sec. 13. 2 1 that galaxies are supported against gravity 
by the fermionic WDM quantum pressure. 

It is highly remarkably that in the context of fermionic WDM the simple static quantum 
description provided by Thomas-Fermi is able to reproduce such broad variety of galaxies. 



Fermions 'at thermal equilibrium ' 
Fermions out of thermal equilibrium 

Observed values taking m = 1 keV 
Observed values taking m = 2 keV 
Observed values taking m = 5 keV 




1.5 



2.5 



3.5 



5.5 



Figure 1: The ordinary logarithm of the square root of the galaxy phase-space density Iog 10 [/12 yQ(P)/m ] = 

3 , , | 

log 10 [34 /}q/(tt ap] defined in eq.{59) as a function of the ordinary logarithm of the product 

4,1 I 

log 10 {»7, (m/kev)3 [pc p(0)/M o ] 6 1 = log 10 [i?o /3q f/,] in parsecs from the numerical resolution of the Thomas-Fermi 
eqs. )50t for WDM fermions. The red (solid) curve is for thermal fermions and the green (dashed) curve corresponds 
to out of equilibrium sterile neutrinos with the distribution eq. jA.6t . The blue crosses + are the observed values of 

K 1 V5(0)/m 2 from Table|2]for m = 1 keV, the red X are the observed values for m = 2 keV and the light blue stars 
are the observed values for m = 5 keV. Notice that the observed Qi, from the stars' velocity dispersion are in fact upper 
bounds for the DM Qi,. 



To conclude, eqs.(l63l indicate that the galaxy magnitudes (halo radius, galaxy masses and 
velocity dispersion) obtained from the Thomas-Fermi quantum treatment for fermion masses in 
the keV scale are fully consistent with all the observations especially for compact dwarfs (see 
Table Namely, fermionic WDM treated quantum mechanically (as it must be) is able to 
reproduce the observed sizes of the DM cores of galaxies. 
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Figure 2: The ordinary logarithm of the galaxy mass log|o[(M/M Q ) ^jM o /[p(0) pc 3 ] (m/keV) 4 ] = 
log 10 [1.425 1 5 A |v'(^o)l / v^o] as a function of the ordinary logarithm of the product 

4 t I i 

log 10 {i, (»j/kev)3 [pc p(0)/M G ]5| = log| [i?o /3g f/,] in parsecs from the numerical resolution of the Thomas-Fermi 
eqs. <5()t for WDM fermions. The red (solid) curve is for thermal fermions and the green (dashed) curve corresponds 
to out of equilibrium sterile neutrinos with the distribution eq. jA.6t . The blue crosses + are the observed values of 
log 10 [(M/M o ) VM /[p(O) pc 3 ] (m/keV) 4 ] from Table|2]for m = 1 keV, the red X are the observed values for m = 2 
keV and the light blue stars are the observed values for m = 5 keV. 



Appendix A. Halo bounds from the Pauli Principle 

We derive here bounds on the phase space density Q(r) from the Pauli Principle assuming a 
simple factorized form for the phase space distribution function 

mp) = n{f) ¥(p/po), 

where po is the momentum characteristic scale of the DM particles. For fermions at thermal 
equilibrium po coincides with the temperature. Since the DM number density follows integrating 
f(r, p) over p as in eq.dlOb. the function ^(p/ po) must fulfill 

J P {2nhf 2n>& J„ d K ) 
We can thus write the phase space distribution function as 

m p ) = n( r> ) l^.JP.) , (A.l) 
P \PoJ 

where the function 

p 3 

2 7T n? po 
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is normalized by 

x 2 dx if/(x) = 1 . 



The Pauli principle imposes the bound eq.© which combined with eq.( dA.U yields, 

«(F) < 



ft 3 7T 2 lfrip/p ) ■ 

Since the momentum distribution ifr(p/po) is normally a monotonically decreasing function of p, 
the most stringent bound is obtained setting p — 0: 

3 ,/T7 -.3 _3 



Po V27 



m cr 



" ft 3 it <A(0) ~ ft 3 7T? <ft0) ■ (A - 2) 



Therefore, the phase space density eq.([T3l is bounded by 



4 

m 



Q(r)<K — , (A3) 



where the dimensionless quantity K is given by 

V27 
n 2 if/(0) 

The value of K depends on the nature of the momentum distribution 



K^^TT^ (A.4) 



For DM particles in thermal equilibrium po is equal to the temperature T and 

T(JC) = -J-r . (A.5) 
e x + 1 



For sterile neutrinos in the v-MSSM model Il52ll_w hich decouple out of equilibrium. Their 
freezed-out distribution function *F(jc) is given by 18[] 

W(x) = 2t /-y _— . (A.6) 



where r » 0.03 is a coupling constant. This formula is valid for all y[x> t and we take 'P(x) = 1 
for s/x < t. 

For sterile neutrinos in the Dodelson-Widrow model 112 ill we have (approximately) the freezed- 
out distribution function 

f 1 

Y(jc) = — — where f * 0.043 keV . (A.7) 

m e + 1 

We display K in Table Q] for thermal fermions, out of thermal equilibrium fermions in the 
DW and v-MSM models and for the Maxwell-Boltzmann distribution. 

Notice that the quantum bound is more restrictive (K turns to be smaller) in the out of thermal 
equilibrium v-MSM model than in the other cases. This is due to the fact that generically, out of 
equilibrium distributions as the v-MSM have more particles with low momentum iSQjJ. if/(Q) is 
therefore larger than at thermal equilibrium and from eq. (IA.4l ) K is smaller. 

The equilibrium Fermi-Dirac distribution produces the less restrictive bound in Table [TJbe- 
cause, somehow, it already knows about the exclusion principle. 
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